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INTRODUCTION  '  . 

■v  -  \  ' 

"The  propagation  of  ultrasound  in  fiber-reinforced  composites  is  controlled  by  the 
relative  magnitudes  of  characteristic  length  variables.  These  length  variables  are  the 

wavelength  (A),  the  fiber  diameter  (a),  the  thicknesses  of  plies  (h),  and  the  overall  thickness 

of  the  component  (H).  It  may  generally  be  assumed  that  A  is  much  larger  than  the  fiber 

diameter  (a  ~  8  pm).  On  the  other  hand  wavelengths  of  the  order  of  thicknesses  of  plies  (h  - 

100  pm),  but  smaller  than  the  overall  thickness  of  the  component  may  well  occur.  In  that 
range  of  wavelengths  the  propagation  of  ultrasonic  waves  gives  rise  to  interesting  multiple 
reflection  phenomena  which  are  the  topic  of  this  paper. 

In  this  paper  we  consider  the  propagation  of  longitudinal  waves  in  a  cross-ply  fiber- 

reinforced  composite.  Since  A  »  a,  we  will  use  an  effective  modulus  representation  for  the 
mechanical  behavior  of  the  plies.  This  means  that  in  the  direction  normal  to  the  plies 
longitudinal  waves  will  propagate  with  a  wave  speed  which  is  independent  of  the  fiber 
direction.  Hence  the  problem  considered  here  is  a  one  dimensional  problem  of  wave 
propagation  through  layers  with  identical  mechanical  impedances  which  are  separated  by 
interfaces. 

Since  the  interfaces  of  the  plies  and/or  laminates  may  reflect  sound,  a  composite  may 
display  some  of  the  features  that  are  characteristic  of  the  propagation  of  ultrasound  in  periodic 
media,  such  as  passing  and  stopping  bands  in  the  frequency  spectrum.  Strict  structural 
periodicity  is,  however,  an  idealization.  In  reality  there  will  be  random  deviations  from 
periodicity,  which  will  in  fact  destroy  the  passing  bands,  and  will  give  rise  to  attenuation  at  all 
frequencies.  This  attenuation  is  often  referred  to  as  localization.  The  present  paper  analyzes 
such  localization.  ( j 

INTERFACES 

Reflection  and  transmission  of  ultrasonic  wave  motion  by  an  interface  or  an  interphase 
between  two  solid  materials  has  been  considered  for  a  number  of  special  cases.  The  insen  of 
Fig.  1  shows  an  incident  wave:  u  =  exp(ikx),  a  reflected  wave:  u  =  R  exp(-ikx),  and  a 
transmitted  wave:  u  =  T  exp(ikx).  Here  we  will  consider  the  case  that  the  materials  on  both 
sides  of  the  interface  or  interphase  are  identical.  Typical  examples  of  interphases  for  such  a 
configuration  are  a  distribution  of  cracks,  as  considered  by  Sotiropoulos  and  Achenbach  [1], 
or  a  distribution  of  spherical  cavities  as  discussed  by  Achenbach  and  Kitahara  [2].  Another 
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Fig.  1.  Magnitudes  of  reflection  and  transmission  coefficients  versus  kh  according  to  Eqn. 
(la,b) 


typical  example  of  an  interphase  is  a  thin  layer  of  slightly  different  material  properties  which 
may  develop  at  a  plane  of  juncture.  An  interphase  is  often  represented  by  a  layer  of  springs. 
For  the  cases  of  cracks,  cavities  and  springs,  the  reflection  and  transmission  coefficients  (R 
and  T,  respectively)  depend  on  the  frequency.  The  general  form  of  the  dependence  of 
absolute  values  on  the  frequency  is  illustrated  in  Fig.  1. 

The  curves  shown  in  Fig.  1  can  conveniently  be  represented  by  functional  forms  that 
depend  on  a  single  dimensionless  parameter  a : 


iakh 
1  -  iakh 


T  = 


1 

1  -  iakh 


(la,b) 


PROPAGATOR  MATRIX  FORMULATION 

Figure  2  shows  a  layer  of  thickness  H,  containing  n  identical  interfaces  separated  by 
distances  h,  thus 

H  =  (n+l)h  (2) 

The  layer  is  coupled  on  both  sides  to  another  medium,  in  this  case  water.  Waves  propagating 
in  the  positive  and  negative  x-directions  have  amplitudes  A  and  B,  respectively. 

First  we  consider  the  interface  j.  Amplitudes  of  waves  to  the  left  of  this  interface  are 
denoted  by  Aj.j  and  Bj.j,  while  amplitudes  to  the  right  are  denoted  by  Ai  and  Bj.  Let  Aj'i 
exp(ikx)  and  Bj  exp(-ikx)  be  considered  as  incident  waves.  Then  if  R  and  T  are  the  reflection 
and  transmission  coefficients,  respectively,  we  have 

Bj.j  —  RAj.i  +  TBj  ,  Aj  =  TAj.i  +  RBj 


(3a,b) 
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Fig.  2.  A  layer  of  darkness  H,  containing  n  identical  interfaces  separated  by  distances  h,  and 
coupled  on  both  sides  to  water 

This  system  may  be  expressed  in  the  form 
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(4a,b) 


where  S  is  the  scattering  matrix.  Since 
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Eqns.  (4)  and  (5)  may  also  be  expressed  as 
/A; 
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where  W  =  S  P  is  the  pi  ^pagator  matrix. 
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Now  let  a  wave  in  the  water  of  amplitude  A.  be  incident  from  the  left  on  the  layer  of 
thickness  H.  The  amplitude  of  the  reflected  wave  is  denoted  by  B-,  and  the  amplitude  of  the 
transmitted  wave  to  the  right  by  A+.  It  is  then  not  difficult  to  see  that  A+,  A.  and  B.  are 
related  by 


fA+\  (A. 

,  =  S0  P  Wn  Wn-i  •••  Wi  So 


B- 
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where 


S0  = 


2k 


w 


/kw+k  kw-k  \ 

V  kw-k  kw+k  J 


(8b) 


Here  So  is  the  scattering  matrix  of  the  solid- water  interface,  kw  is  the  wavenumber  in  water, 
kw  =  CO/Cw. 
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Equation  (8a)  can  be  solved  for  A+  and  B.  in  terms  of  A,,  the  amplitude  of  the  incident 
wave.  The  reflection  and  transmission  coefficients  of  the  layer  may  then  be  defined  as 


(9a, b) 


For  the  case  that  all  interfaces  are  the  same  and  of  the  form  shown  in  Fig.  1,  results  are 

displayed  in  Fig.  3  for  the  case  a  =  0.1  and  n  =  21.  The  peaks  at  kwh  =  0.11  correspond  to 
resonances  for  the  overall  thickness  of  the  layer.  The  drop  of  the  transmission  coefficient  to 
zero  and  the  corresponding  increase  of  the  reflection  coefficient  to  one  occurs  at  kwh  =  4.4, 
which  is  just  the  frequency  at  which  the  stopping  bands  start  for  wave  propagation  in  the 
unbounded  periodic  medium. 


Fig.  3.  Magnitudes  of  reflection  and  transmission  coefficients  of  the  layer  versus  kh,  for  a  = 
0.1,  n=21 

INTERFACES  WITH  STATISTICAL  PROPERTIES 

For  composite  materials  periodicity  is  an  idealization,  since  the  real  ply-structure  will 
always  have  some  random  features.  Within  the  context  of  this  paper's  approach,  statistical 

properties  of  an  interphase  can  be  taken  into  account  by  considering  the  parameter  a  in  Eqn. 
(1)  as  a  random  variable.  The  propagation  in  a  basic  cell  can  still  be  defined  by  a  propagator 
matrix,  but  now  this  matrix  has  random  coefficients.  Following  the  approach  developed  for 
the  deterministic  system,  see  Eqn.  (8a, b),  we  have  to  consider  the  product  of  a  finite  number 
of  random  matrices.  Unfortunately  it  is  very  difficult  to  make  rigorously  valid  statements 
about  such  a  product.  It  is,  however,  possible  to  draw  rigorously  valid  conclusions  when  the 
number  of  matrices  becomes  very  large.  This  can  be  done  by  the  use  of  a  theorem  due  to 
Furstenberg  [3]. 
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When  the  number  of  interfaces  is  large,  the  asymptotic  behavior  of  the  total 
transmission  coefficient  T /  has  a  close  relationship  with  the  magnitude  of  eigenvalues  of  the 

propagator  matrix  W.  Let  exp(X)  be  the  eigenvalue  of  W,  then  lexp(X)l  =  exp[Re(X)].  If 
Rc(X)  =  0,  as  in  passing  bands  of  the  periodic  structure,  lexp(X)l  =  1,  and 

t/  =  ^7~°(1)  (10) 

If,  on  the  other  hand,  Re(X)  *  0,  as  in  stopping  bands  of  the  periodic  structure,  lexp(X)l  *  1, 
and 


T/  =  ^~o(l)  (11) 

Furstenberg’s  theorem  states  that  for  the  mean  value  of  the  eigenvalues,  exp(V),  of  the 
random  matrix  W,  Re(X)  will  never  be  equal  to  zero,  which  implies  that  for  a  solid  with 
random  interfaces,  at  all  frequencies  we  will  have 

T/~o(l)  (12) 

Since  in  our  problem,  conservation  of  rates  of  energies  holds,  the  information  on  T / 
subsequently  supplies  the  information  on  R/. 

Filrstenberg's  theorem  may  be  stated  as  follows:  For  random  matrices  Wj,  i  =  1,2, 

•  •  •,  n,  there  exists  a  y  >  0,  such  that  for  each  Zo  *  0, 

lim  ilnlWn--W2Wizol  =  y  (13) 
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with  probability  of  unity,  and 
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(15) 


Here  a  and  0  are  random  variables  in  W  and  zo,  respectively,  and  p(a)  and  v(0)  are 

corresponding  probability  density  functions.  The  constant  y  is  usually  called  the  localization 
constant,  while  /  =  1/y  is  called  the  localization  length  [4]. 

Suppose  that  the  phase  of  zo  is  evenly  distributed  over  the  interval  [0,  tc],  and  the 
interface  parameter  a  obeys  a  Gaussian  distribution.  Then,  assuming  that  a  deviates  slightly 
from  its  mean  value,  it  can  be  shown  by  evaluating  Eqn.  (15)  that 


y«  Oa  co2 
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where  aa  is  the  deviation  of  the  interface  parameter  a.  Let 
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then  according  to  Fiirstenberg  theorem,  ay  ~  exp(ny).  By  solving  Eqn.  (8)  for  T/  and  R /,  and 
using  Eqn.  (18),  we  obtain 


T i  = 


A± 

A- 


—  -  exp(-ny) 


(19) 


(It  can  be  proven  that  det(ay)  =  1). 


Fig.  4.  Localization  constant  y  versus  kh  for  a  =  0.1,  ca  =  0.3 


Figures  4  and  5  show  the  localization  constant  y  and  the  localization  length  /  as  functions 
the  dimensionless  wave  number  kh.  Results  based  on  Fiirstenberg  theorem  aie  compared 
with  those  obtained  by  the  use  of  Monte  Carlo  simulation  and  for  the  stnctly  periodic 
structure.  These  results  have  been  computed  for 


&  =  0.1  and  On  =  0.3 


(20) 
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Fig.  7.  Magnitude  of  transmission  coefficient  versus  kh  for  a  layer  of  200  interfaces,  a  = 

0.1,  aa  =  0.3 

where  the  bar  denotes  the  mean.  It  is  important  to  note  that  y  is  not  equal  to  zero  in  the  range 
of  kh  which  would  correspond  to  the  passing  bands  of  the  strictly  periodic  composite.  The 
results  presented  in  Fig.  6  and  Fig.  7  show  the  direct  matrix  multiplication  for  T/  and  R /, 
where  the  number  of  interfaces  is  200.  It  can  be  easily  seen  that  the  transmission  coefficient 
T/  decreases  with  increasing  frequency  and  correspondingly,  the  reflection  coefficient  R/ 
increases  with  increasing  frequency. 


CONCLUSION 

It  has  been  shown  in  this  paper  that  the  propagator  matrix  approach  can  be  used  to 
study  one  dimensional  propagation  of  ultrasound  in  composites.  Random  variations  of  the 
interface  properties  give  rise  to  the  localization  phenomenon.  As  a  consequence,  ultrasound 
will  be  confined  to  a  local  region  near  the  insonified  area.  It  should  be  expected  that 
localization  will  adversely  affect  ultrasonic  testing  of  thick  composites. 
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In  this  paper  we  consider  the  propagation  of  longitudinal  waves  in  a  cross-ply 
r iber-re inf orced  composite.  Since  A  >  >  a,  we  will  use  an  effective  modulus 
representation  for  the  mechanical  behavior  of  the  plies.  This  means  that  in  the 
direction  normal  to  the  plies  longitudinal  waves  will  propagate  with  a  wave  speed 
which  is  independent  of  the  fiber  direction.  Hence  the  problem  considered  here  is 
a  one  dimensional  problem  of  wave  propagation  through  layers  with  identical 
mechanical  impedances  which  are  separated  by  interfaces. 

Since  the  interfaces  of  the  plies  and/or  laminates  may  reflect  sound,  a 
composite  may  display  some  of  the  features  that  are  characteristic  of  the 
propagation  of  ultrasound  in  periodic  media,  such  as  passing  and  stopping  bands 
in  the  frequency  spectrum.  Strict  structural  periodicity  is,  however,  an 
idealization.  In  reality  there  will  be  random  deviations  from  periodicity, 
which  will  in  fact  destroy  the  passing  bands,  and  will  gi,re  rise  to  attenuation 
at  all  frequencies.  This  attenuation  is  often  referred  to  as  localization.  The 
present  paper  analyzes  such  localizat ion . 


